LISA data analysis I: Doppler demodulation 
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The orbital motion of the Laser Interferometer Space Antenna (LISA) produces amplitude, phase 
and frequency modulation of a gravitational wave signal. The modulations have the effect of spread- 
ing a monochromatic gravitational wave signal across a range of frequencies. The modulations en- 
code useful information about the source location and orientation, but they also have the deleterious 
affect of spreading a signal across a wide bandwidth, thereby reducing the strength of the signal 
relative to the instrument noise. We describe a simple method for removing the dominant, Doppler, 
component of the signal modulation. The demodulation reassembles the power from a monochro- 
matic source into a narrow spike, and provides a quick way to determine the sky locations and 
frequencies of the brightest gravitational wave sources. 



I. INTRODUCTION 

This paper introduces a "quick and dirty" method for 
making a first pass at the LISA data analysis problem. 
LISA will be sensitive to gravitational wave sources in 
the frequency range 1CP 4 to 1CP 1 Hz, which nicely com- 
plements the frequencies range covered by the ground 
based detectors (10 1 to 10 3 Hz). The data analysis chal- 
lenges posed by LISA are very different from those en- 
countered with the ground based detectors. Unlike the 
situation faced by the ground based observatories, most 
of the sources that LISA hopes to detect have well un- 
derstood gravitational waveforms - indeed many of them 
will be monochromatic and unchanging over the life of 
the experiment. The difficulty comes when we include 
LISA's orbital motion and the modulations this intro- 
duces into the signal. Since LISA is designed to orbit at 
I AU, the modulations introduce sidebands at multiples 
of the orbital frequency of f m = I/year. 

The approach we investigate here works by correct- 
ing for the dominant, Doppler modulation. This has the 
effect of re-assembling most of the Fourier power of a 
monochromatic source into a single spike. Since the cor- 
rection depends on where the source is located on the 
sky, we are able to locate a source and its frequency. 

The Doppler demodulation approach is familiar to ra- 
dio astronomers, and has been discussed in relation to 
gravitational wave astronomy H. Doppler demodulation 
forms an integral part of ground based gravitational wave 
searches for continuous gravitational wave sources such as 
pulsars (2[ [|. However, the unique orbital motion of the 
LISA observatory, and the different frequency range that 
it covers, demand a separate treatment of the Doppler 
demodulation approach to LISA data analysis. 



II. SIGNAL MODULATION 

LISA will output a time series describing the strain in 
the detector as a function of time, h(t). The strain will 



be a combination of the signal, s(t), and noise, n(t), in 
the detector: h(t) = s(t) + n{t). Because of LISA's or- 
bital motion, a source that is monochromatic at the Sun's 
barycenter will be spread over a range of frequencies. A 
monochromatic gravitational wave with frequency / and 
amplitudes h+ and h x in the plus and cross polarizations 
will produce the response m 



where 



s(t) = A(t) cos ^>{t) (I) 



#(t) = [2wft + <p + <f> D {t) + <f>p(t)} • (2) 



The amplitude modulation A{t), frequency modulation 
4>o{t) and phase modulation 4>p(t) are given by 
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A(t) = [(h + F+(t)) 2 + (h x F*(t)) 



R 

b D (t) = 2irf—sm9 s cos(2Trf m t-(j> s ) 
c 



bp(t) = — arctan 



h x F*(t) 
h+F+(t) 



(3) 



(4) 



(5) 



Here F + (t) and F x (t) are the detector antenna patterns 
in barycentric coordinates Each of the modulations 
is periodic in the orbital frequency f m = 1/ycar. The 
angles 9 S and 4> s give the sky location of the source, ifo 
is the phase of the wave and R — 1 AU is the Earth-Sun 
distance. The frequency modulation, which is caused by 
the time dependent Doppler shifting of the gravitational 
wave in the rest frame of the detector, is the dominant 
effect for frequencies above 5 x ICU 4 Hz||. Ignoring the 
amplitude and phase modulation for now, pure Doppler 
modulation takes the form 

s(t) = Acos [2irft + {3 cos(27r f m t — 4> s ) + tpo] 

= $(a J2 J n (/3)e 2 ^ {f+f '" n)t e lVo e m ^ /2 -^AQ) 
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Here J n is a Bessel function of the first kind of order n 
and (3 is the modulation index 



R 

13 = 2nf— sin( 

c 



(7) 



The bandwidth of the signal - defined to be the frequency 
interval that contains 98% of the total power - is given 
by 



5 = 2(1 + /?)/„ 



(8) 



Sources in the equatorial plane have bandwidths ranging 
fromB = 2.6xl(T 7 Hz at / = 1CT 3 Hz to B = 2.1 xl(T 6 
Hz at / = f(T 2 Hz. 



III. DEMODULATION 

If LISA were at rest with respect to the sky, each 
monochromatic source would produce a spike in the 
power spectrum of s(t). While this would make data 
analysis very easy, it would also severely limit the sci- 
ence that could be done as it would be impossible to 
determine the source location or orientation. Since LISA 
will move with respect to the sky, each source will have 
its own unique modulation pattern, and this pattern can 
be used to fix its location and orientation. On the other 
hand, the modulation makes the data analysis more com- 
plicated as the Fourier power of each source is spread 
over a wide bandwidth B. One approach to the data 
analysis problem is to demodulate the signal, thereby 
re-assembling all the power of a given source at one fre- 
quency. Since each source has a unique modulation, each 
source will also require a unique demodulation. Demod- 
ulating a particular source is not difficult if one happens 
to know its frequency, location, orientation and orbital 
phase. One could imagine searching through this param- 
eter space and looking for spikes in the power spectrum 
corresponding to sources that have been properly demod- 
ulated. A more practical approach is to focus on the 
Dopplcr modulation as it causes the largest spreading of 
the Fourier power. 

Consider the Doppler modulated phase 



*(t) = 2tt/ 



t H sin 6 S cos(27r/ m i — <p s 

c 



+ ¥>o • (9) 



We seek a new time coordinate t' in which this phase is 
stationary: 



d<P d$ dt 

dF ~ ~dtdF 



(10) 



1 sin# s sin(27r/ m i — <j> s ) ) dt . (11) 



Thus, 



Working to first order in v/c we have 
R 

t = t' sin0 s cos(27r/ m t-& i ) (12) 

c 



Taking the data stream from the detector over a one 
year period and performing a fast Fourier transform al- 
lows us to write 



s(t) = ^2 a n& 



27rif rn nt 



(13) 



Performing the coordinate transformation ( |12[ ) we arrive 
at the new Fourier expansion 



27rif m kt 



<0 = 5>*< 



where the Fourier coefficients c k are given by 



(14) 



R 



Ck 



= y j a n J n ^ k (2nnf m -sm0 s )e^ n ~ k ^-^ 2 ^ . (15) 



Since the modulation has a limited bandwidth, an excel- 
lent approximation to c k is given by 



k+l 



Ck — a n J n - k (ci)e lin ~ 



fc)0„-7r/2) 



=k-l 



where 



and 



R , 

a = ZTTKjm — sin ( 
c 



I = 1 + la] 



(16) 



(17) 



(18) 



The square brackets imply that we have taken the in- 
teger part of a. Consider the pure Doppler modulation 
described in (^|) for the simple case where / = qf m , and 
q is an integer (the non-integer case is more involved, but 
the idea is the same). The source has Fourier components 

a„ = Ae^^^e^'W 2 -*-' , (19) 

so that 

c k -ie^e'^)*-/^ J n ^ q (0)J n ^ k (a) . (20) 

n 

Using the Neumann addition theorem, 

Jn- q (P)Jn-k{a) = J k - q {a - /?), (21) 



and the fact that a sa /3, we find 

Ck~Ae^S kq . 



(22) 



In other words, the Doppler demodulation procedure re- 
assembles all the power at the barycenter frequency / = 
qfm- The demodulation is less effective when applied to 
a real LISA source because it does not correct for the 
amplitude and phase modulations. 

Our Fourier space approach is very efficient if one is 
interested in a limited frequency range. It does not of- 
fer any great saving over a direct implementation in the 
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FIG. 1: Power spectra, P(f), before (dashed line) and after 
(solid line) Doppler demodulation. 





/ (Hz) 


6s 


4> s 


A 


4.999904 x 10" 3 


2.867424 


0.258971 


B 


4.999729 x 10" 3 


0.661738 


3.323842 



TABLE I: Description of Sources A and B 



time domain if one wants to consider all frequencies at 
once. HellingsQ has recently implemented the Doppler 
demodulation procedure in the time domain and found 
similar results to ours. 

To illustrate how the Doppler demodulation works we 
consider monochromatic, circular Newtonian binaries as 
described in Ref. Figure 1 shows the power spec- 
trum of s(t) before and after Doppler demodulation for 
a source with / = 5.000167 x 10~ 3 Hz, 9 S = 2.385366 
and 4> s — 4.462868. We see that most of the power is col- 
lected into a spike of width ~ 3/ m about the barycenter 
frequency. 

Next we consider an example where there are two 
sources that are nearby in frequency, but nevertheless do 
not have overlapping bandwidths. The sources have / = 
5.001419 x 10~ 3 Hz, (0 S)( £ S ) = (2.056685,2.212044) and 
/ = 4.9988 x 10~ 3 Hz, (6» S ,0 S ) = (1.280671,1.446971). 
The higher frequency source has an amplitude 1.103 
times larger than the lower frequency source. The effect 
of demodulating each source is shown in Figure 2. 

The performance of the demodulation procedure is 
considerably less impressive when the two sources have 
overlapping bandwidths. Consider sources A and B de- 
scribed in Table 1. The amplitude of Source B is 1.59 
times larger than the amplitude of Source A. Figure 3 
shows what happens when we demodulate Source A and 
Source B in turn. Because the two sources overlap, their 
signals interfere with each other, and the demodulation 
is only able to detect the stronger of the two sources. 
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FIG. 2: Power spectra before (dashed line) and after (solid 
line) Doppler demodulation. The upper graph shows the de- 
modulation when the detector is made stationary with respect 
to the higher frequency source, while the lower graph shows 
the same for the lower frequency source. 



A. Angular resolution 

The Doppler demodulation technique can be used to 
find the sky location of a source by successively demod- 
ulating each point on the sky. In practice we use the 
HEALPIXj?) hierarchical, equal area pixelization scheme 
to provide a finite number of sky locations to demodu- 
late. We search for the maximum power contained in 
three adjacent frequency bins, and record this value for 
each sky pixel. Rather than find the maxima across all 
frequencies, we produce separate sky maps for frequency 
intervals of width 47r/^/ m . There is little point in trying 
to use smaller frequency intervals due to the interference 
effect discussed earlier. 

We begin by considering an isolated source (i.e. no 
other sources with overlapping bandwidths) with / = 
5.00015 x 10~ 3 Hz and (9 S ,4> S ) = (1.851478,4.934289). 
Using sky pixels of angular size ~ 0.92° we arrive at 
the graph shown in Figure 4. The demodulation code 
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FIG. 3: Power spectra before (dashed line) and after (solid 
line) Doppler demodulation. The upper graph shows the de- 
modulation when the detector is made stationary with respect 
to Source A, while the lower graph shows the same for Source 
B. 



producedthe best fit values of /_= 5.00014 x 10 3 Hz 
and {6 s ,cf> s ) = (1.253,5.04) and (0 S ,0 S ) = (1.889,5.04). 
The degenerate fit for the sky location illustrates one 
of the drawbacks of the Doppler demodulation method 
- it is unable to distinguish between sources above or 
below the equator. The errors in the source location, 
(A9 S — 2.1°, A<fi s — 6.3°), are larger than the pixel size, 
and can be attributed to the amplitude and phase mod- 
ulations which we have not corrected for. It should be 
noted that the error in the source location is not caused 
by instrument noise since we are working in the large 
signal-to-noise limit (n(t) = 0). Indeed, instrument noise 
has little effect on the demodulation procedure since the 
noise is incoherent. The noise gets moved about in fre- 
quency space, but there is little power accumulation at 
any one frequency. When noise was added to the previous 
example at a signal to noise ratio of 5, the frequency de- 
termination and 4> s determination were unaffected, while 
the error in 9 S grew to A9 S = 3.4°. 



FIG. 



4: Determining the sky location of a source with 
(9 S ,4> S ) = (1.851478,4.934289). The sky map uses the Moll- 
weide projection in Ecliptic coordinates with (n/2,0) at the 
center of the map. 



FIG. 5: Determining the sky location of source A without 
source B. 



Finally, we illustrate how the source interference phe- 
nomena discussed earlier limits our ability to locate 
sources that overlap in frequency. Figures 5 and 6 show 
how the demodulation procedure is able to locate Sources 
A and B when one of the sources is turned off. Figure 
7 shows what happens when both sources are present. 
We see that source B shows up clearly, while source A is 
washed out. 



IV. DISCUSSION 

The Doppler demodulation procedure provides a quick 
way of finding the frequencies and sky locations of the 
brightest sources detected by LISA. The method has a 
number of limitations, the most serious being its inabil- 
ity to locate more than one source per bandwidth, and 
its inability to determine if a source is in the northern 
or southern hemisphere. Despite these limitations, the 
Doppler demodulation procedure will be a useful tool in 
the LISA data analysis arsenal. 



5 





FIG. 7: Determining the sky location of sources A and B 
together. 
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